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SUMMARY 


The design of nonlinear controllers has relied on the use of detailed aerodynamic and engine models that must 
be associated with the control law in the flight system implementation. Many of these controllers have been applied 
to vehicle flightpath control problems and have attempted to combine both inner- and outer-loop control functions 
in a single controller. In this paper, a new approach to the nonlinear trajectory control problem is presented. This 
approach uses linearizing transformations with measurement feedback to eliminate the need for detailed aircraft 
models in outer-loop control applications. By applying this approach and separating the inner-loop and outer-loop 
functions two things are achieved: (1) the need for incorporating detailed aerodynamic models in the controller is 
obviated and (2) the controller is more easily incorporated into existing aircraft flight-control systems. This paper 
also discusses an implementation of the controller described here. This controller is tested on a six degree-of- freedom 
F-15 simulation and in flight on an F-15 aircraft. Simulation data are presented which validates this approach over 
a large portion of the F-15 flight envelope. Proof of this concept is provided by flight-test data that closely matches 
simulation results. Flight-test data are also presented. 

INTRODUCTION 

The problem of aircraft trajectory control is generally thought of as an “outer-loop” problem, that is, the air- 
craft is assumed to have a stabilizing, “inner-loop” command augmentation system through which the pilot or au- 
tomatic trajectory control system interacts with the vehicle control effectors (such as the elevators, ailerons, and 
rudder). In fact, for a great number of trajectory control problems, it is the pilot who provides this outer-loop 
control by translating a mental image of the trajectory into stick, pedal, and throttle commands using feedback in- 
formation provided by the crew station instruments. Outer-loop control may also be provided by an autopilot that 
performs simple tasks like maintaining altitude and velocity or one that performs more complex tasks such as auto- 
matic landing. 

Inner-loop control generally refers to control of the aircraft about its center of mass; outer-loop control refers to 
the control of the trajectory of the center of mass. The former is generally called “the control problem” while the 
latter is often referred to as “the guidance problem.” 

This distinction between inner-loop and outer-loop functions is important for two reasons: (1) the frequency 
response requirements of these functions are significantly different and (2) the output commands from the two types 
of controllers differ significantly. 

The frequency response requirements for inner-loop control are generally on the order of 1 to 30 rad/sec, because 
of the dynamic characteristics of high performance aircraft. The frequency response requirements for outer-loop 
control are on the order of 0. 1 to 3 rad/sec, because of the trajectory tracking requirements. The outputs from an inner- 
loop control system are essentially surface position commands; the outputs from an outer-loop controller are similar 
to stick, pedal, and throttle position commands. The two types of controllers are usually separate designs and dis- 
tinct subsystems. 

The fundamental objective of trajectory control research at the Dryden Flight Research Facility is to define a 
generic approach for designing control laws that will command an aircraft to fly precise flight-test maneuvers to ob- 
tain accurate, repeatable flight data. The intent is to use modern control theory design techniques to lay out a method 
of designing flight-test trajectory controllers for various aircraft at Dryden to improve data collection performance. 
The approach should also be applicable to other aircraft for the design of items such as autopilots and automatic 
navigation systems. Implicit in this objective is the requirement to develop control laws using different techniques 
for a candidate aircraft and subsequently to evaluate the performance and assess the limitations of these control laws. 
The candidate aircraft is the NASA F-15 highly integrated digital electronic control (HIDEC) aircraft [1], 
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Over the past several years, techniques such as eigenstructure assignment, linear optimal control, and others 
were applied in the attempt to find a design approach. However, the most promising technique involves the use of 
nonlinear inverse transformation theory. Inverse transformations relate the input of the linearized aircraft system to 
the required aircraft command augmentation system (CAS) nonlinear input. The application of this theory to flight- 
test trajectory control has been treated by Menon et al [2,3]. This application uses the equations of motion and inverse 
aerodynamic and engine models to provide inputs to the inner-loop controller of the aircraft. The methodology for 
the use of nonlinear inverse transformations and inverse models presented in this paper is based primarily on the 
work of Meyer and Cicolani [4,5], Meyer [6], and Meyer et al [7] whose approach is presented in some detail. The 
perceived limitations of the work by previous authors are discussed. 

This paper presents the theory for and an application of an approach that provides a new solution to the flight-test 
trajectory control problem. The elegance of this solution is that the control laws are obtained from flat, nonrotating 
Earth, equations of motion with accelerometer outputs providing force data instead of an aerodynamic model. This 
makes the controller easy to apply to different aircraft since only a simple inverse engine model is required. The use 
of nonlinear inverse transformations presents a linear system to the designer for which simple proportional-integral, 
or proportional-integral-derivative control laws using a single set of gains need be designed. 

A flight-test trajectory controller (FTTC) was designed using this new approach and was tested on a real-time, 
six-degree-of-freedom, full envelope, nonlinear simulation of an F-15 aircraft. Various maneuvers were simulated 
over a Mach number range of 0.60 to 1.20, over an altitude range of 10,000 ft to 40,000 ft. Simulation results justified 
flight test of the FTTC with the F-15 HIDEC aircraft. Flight test covered a subset of the maneuvers tested in the 
simulation. Both simulation and flight-test results are presented. 

NOMENCLATURE 


Acronyms 


AWJSRA 

augmentor wing jet short takeoff and landing research aircraft 

CAS 

command augmentation system 

FTTC 

flight-test trajectory control 

FTTG 

flight-test trajectory guidance 

HIDEC 

highly integrated digital electronic control 

HiMAT 

highly maneuverable aircraft technology 

PLA 

power lever angle, deg 

RAV 

remotely augmented vehicle 

STO 

short takeoff and landing 

TAFCOS 

total automatic flight-control system 

VATO 

vertical attitude takeoff and landing 

WATR 

Western Aeronautical Test Range 


Symbols 

A 

angles; plant matrix 

A o 

servo plant matrix 

a 

acoustic velocity, ft/sec 
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a z 
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PLA com 
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Pcom 
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Qcom 
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r* 

SR 

S 
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u 

u 

V 
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body axis normal accelerometer output, g 

x-axis body accelerometer output, g 

y-axis body accelerometer output, g 

z-axis body accelerometer output, g 

control matrix 

servo control matrix 

aircraft total drag, lb 

error signal vector 

error signal 

aircraft nonlinear functions 
aircraft nonlinear control functions 
gravitational acceleration, ft/sec 2 
aircraft measurement functions 
altitude, ft 
altitude rate, ft/sec 
vertical acceleration, ft/sec 2 
gain matrix 

total aerodynamic lift, lb 
Mach number 
mass, slugs 

commanded power lever angle, deg 
body axis roll rate, dcg/sec; power state 
power rate, ft/sec 

commanded body axis roll rate, deg/sec 
body axis pitch rate, deg/sec 
commanded body axis pitch rate, deg/sec 
body axis yaw rate, deg/sec; position state 
reference trajectory 
set of real numbers 

transformation from body-axis rotational rates to inertial angle rotational rates 

linearizing transformation 

total thrust, lb 

nonlinear input vector 

aircraft x-body axis velocity, ft/sec 

aircraft x-body axis kinematic acceleration, ft/sec 2 

total velocity, ft/sec 

total kinematic acceleration, ft/sec 2 


linear input vector 

aircraft y-body axis velocity, ft/sec 

aircraft y-body axis kinematic acceleration, ft/sec 2 

inverse linearizing transformation 

noise 

aircraft z-body axis velocity, ft/sec, noise 
aircraft z-body axis kinematic acceleration, ft/sec 2 
thrust in the aircraft x-body axis, lb 
state vector 

state time derivative vector 

state vector of the apparent linear plant 

north inertial position, ft 

aircraft total sideforce, lb 

measurement vector 

east inertial position, ft 

thrust in the aircraft z-body axis, lb 

vertical inertial position, ft 

angle of attack, deg 

angle of attack rate, deg/sec 

angle of sideslip, deg 

wind reference frame Euler pitch angle, deg 
commanded incremental load factor, g 
engine orientation angle, deg 
throttle, deg 

commanded throttle, deg 

body reference frame Euler pitch angle, deg 

body reference frame Euler pitch angle rate, deg/sec 

body reference frame Euler roll angle, deg 

body reference frame Euler roll angle rate, deg/sec 

rotational rates, deg/sec 

rotational accelerations, deg/sec 2 


Subscripts 


aerodynamic 

body axis normal accelerometer output 
body axis 

body axis relative to inertial system 



D 

derivative 

h 

altitude 

I 

integral 

k 

kinematic 

M 

Mach number 

m 

measurement feedback system 

P 

proportional 

P 

power-configuration model 

P 

nonlinear plant system; power 

Q 

body axis pitch rate 

r 

position 

ref 

reference value 

u 

function of the input 

V 

total velocity 

V 

estimate or model output 

X 

function of the state 

y 

output 

0 

initial conditions, or zero altitude 

a 

angle of attack 

4 > 

body reference frame Euler roll angle 


BACKGROUND 

The motivation for using direct nonlinear flight-control synthesis techniques for designing aircraft autopilots 
and guidance systems is that these techniques are not limited by their reliance on linear perturbation models of 
the aircraft. While linear design techniques easily handle linear or nearly linear systems, as system nonlinearities 
become more pronounced, the effort required to design guidance laws increases rapidly. Although linear techniques 
are still capable of addressing such problems, a new technique is desired that would handle the nonlinearities simply 
and more directly. 

The flight-test trajectory control problem [8,9,10,11] is a particular type of guidance problem that faces these 
issues. Flight-test maneuvers often cover ranges of the flight envelope, encountering aerodynamic and thrust non- 
linearities. In some cases, the maneuvers specified are themselves nonlinear trajectories such as windup turns and 
pushover— pullups. These effects increase the complexity required from a linear design technique. 

Flight-Test Trajectory Control 


In flight-test trajectory control, the requirement is to provide precise, repeatable control of an aircraft during ma- 
neuvers used to gather aerodynamic, structural, propulsion, and performance data. For conventional flight regimes 
these maneuvers are fairly straightforward and, in general, easily modeled. Level-accelerations-decelerations, 
pushover— pullups, and windup turns are typical of these maneuvers. The original motivation in developing pilot- 
aiding displays for the open-loop problem [9] and closed-loop controllers for automatic maneuvering [8,10,1 1 ] was 
to increase the quality of flight data and to ensure repeatable, consistent maneuvers. 
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Numerous approaches to the flight-test trajectory problem were taken including cut-and-try techniques on a non- 
linear simulation [9], linear optimal control [12,13,14], optimal cooperative control [15], eigenstructure assignment 
[13,16,17], suboptimal output feedback [12], and classical control theory [11]. Although successful open-loop and 
closed-loop systems were developed and used in flight, each of these approaches had limitations. The open-loop 
flight-test trajectory guidance (FTTG) was used on a number of aircraft programs [9] and was developed using a 
cut-and-try approach on piloted, nonlinear simulations. The success of the FTTG suggested that even better results 
could be obtained using a closed-loop controller, and the flight-test maneuver autopilot (FTMAP) was developed 
for the highly maneuverable aircraft technology (HiMAT) vehicle using classical control methods [8,10,11], How- 
ever, while there were successful applications of solutions to the flight-test trajectory control problem, none of the 
design techniques employed in the development of these systems were adequate in terms of providing a concise 
methodology for developing outer-loop controllers. 

Except with the FTTG, all of the problems of transforming from a linear to a nonlinear system were encountered 
and an extensive cut-and-try session was needed to adjust the controllers when they were implemented and tuned 
using a nonlinear simulation. Furthermore, because the most difficult flight-test maneuvers have significant coupling 
between axes, the tuning problem was more than simply the iterative adjustment of parameters within a single control 
loop; simultaneous adjustment of parameters in multiple control loops was required. The use of nonlinear control 
[2,3] simplifies many of these problems, although some cut-and-try tuning is still required. 

Deriving linear perturbation models of aerodynamics and propulsion systems [ 1 8,19] as well as inner-loop, flight- 
control systems is a difficulty facing the engineer using design methods based on linear perturbations. Perhaps an 
even more formidable problem that arises in using linear perturbation methods for FTTCs is the requirement for a 
large number of design cases and, hence, perturbation models. In addition to the requirements imposed because aero- 
dynamic forces and moments change as a function of flight conditions, varying degrees of cross-axis coupling often 
must be accounted for at each flight condition. These additional requirements add significantly to an already complex 
design problem. 


Nonlinear Control 

Meyer and Cicolani [4] describe a formal structure for advanced flight-control systems that provides the basis 
for the results described in this paper. This formal structure may now be taken as preliminary results that motivated 
a large body of work in nonlinear controls and inverse model theory. This formal structure included features that 
were assumed in later applications: 

1. A block triangular system was assumed, 

2. inner- and outer-loop control functions were not explicitly separated, 

3. a time-scale separation (hierarchy) was assumed (either explicitly or implicitly), 

4. inverse aerodynamic and thrust models were used to provide an “open-loop, a priori" feedforward control, and 

5. the feedforward control included an inverse model of the aircraft and control system dynamics. 

This formal structure, the total automatic flight-control system (TAFCOS), has been applied in detailed nonlinear 
simulations of the de Havilland C-8A Buffalo augmentor wing jet short takeoff and landing (STO) research aircraft 
(AWJSRA) [20,21], and a vertical attitude takeoff and landing (VATO) aircraft [22,23,24], The TAFCOS has also 
been used in flight tests on a DHC-6 Twin Otter aircraft [25,26], the AWJSRA [5], and a UH-1H helicopter [7,27], 
Excellent results were achieved in all of these studies and the claim made by Meyer et al [27] that “the method is ef- 
fective for a large class of dynamic systems that require multi axis control and that have highly coupled nonlinearities, 
redundant controls, and complex multidimensional operational envelopes” can be taken as thoroughly validated. 
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Parallel to the cited work on TAFCOS, Meyer, Hunt, and Su [6,7,28,29,30,31,32,33] have pursued research in 
the underlying theory of linearizing transformations by model inversion. Necessary and sufficient conditions for the 
existence of transformations from nonlinear to linear systems were established by Su [32] and Hunt et al [30]. These 
transformations were shown to be global results by Hunt et al [30], Meyer [6] and Meyer et al [7] have described a 
general structure (fig. (1)) for automatic control based on these results. This structure contains transformations (the 
WT-map) which when taken with the aircraft, provide a linear system that is modeled in the ‘MODEL SERVO. The 
reference trajectory, r\ is fed into the ‘MODEL LAW’ which determines the required input for the linear system to 
follow that trajectory. The regulator then determines the input required for the W-map using the linear input, v 0 , the 
linear states, y 0 , and the apparent linear system state feedback, y. 



Model servo 


Regulator WT-map 


Aircraft 

900350 


Figure 1. Meyer’s formal structure for advanced flight-control systems. 


In this structure, a general nonlinear system equation of the form 

x(t) =f[x(t),u(t)] O) 

with the state x e 5R n , the control u e SR™, and f : 3R n x 5R m -> $R n , an analytic function is assumed. Because it is 
a necessary condition for such a system to be transformable to a controllable linear system [27], it is assumed that 
equation (1) can be reformulated as an equivalent system in terms of new state and control coordinates x and u so 
that the system is in the form 

m 

x(t) = f [x(t)] + 52 u i(t)gi[x(t)] (2) 

«=i 


with 


f CO) = 0 


(3) 


This nonlinear system is transformed using a mapping T : 5R n — ♦ SR* x 5R* , called the T-map, such that 


y 

v 


r(x,u) 


(4) 
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with y 6 5R , v £ 3^. where y and v are consistent with the state and control vectors of a time-invariant, linear 
system of the form 

yo(0 = Aoy 0 (t) + Sov 0 (0 (5) 

This linear system model is used to determine reference inputs, vo, into a form compatible with the linear con- 
trol law component of the linear controller (“REGULATOR” in fig. (1)). The linear system in equation (5) is 
sometimes referred to as the reference model” and the elements of v in the linear, transformed system are known 
as “pseudocontrols.” 

An inverse T-map (sometimes referred to as a W-map) maps the pseudocontrols of the linear system, v, into the 
nonlinear system inputs, u, with T~ ] = W : K* x R l -» $R n such that 

u = W'Xy.v) (6) 

Also shown in figure 1 is a linear regulator that provides feedback control for the linear system with 

<5v = Ke y ( 7 ) 

v = v 0 + 6v (S) 

ey = y - yo (9) 

The linear reference model, shown as the “MODEL SERVO” in this figure, provides feedforward control for the total 
nonlinear control system. The nonlinear aircraft model to be controlled (“AIRCRAFT’ block in fig. (1)) is presented 
by Meyer [6] and shown in figure 2. In this model, the aircraft is broken into three blocks: power/configuration, 
rotation, and translation. Basically, Meyer’s aircraft model is a full state, six degree-of-freedom model, but has no 
inner-loop control system associated with it. 

A simple example of a nonlinear controller designed to fly a level acceleration test maneuver is derived in 
appendix A. 
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Figure 2. Meyer’s generic aircraft model. 

Application of Nonlinear Control to the Flight-Test Trajectory Control Problem 

Menon et al [2] report the development of a nonlinear controller for the F-15 aircraft flight-test trajectory control 
problem. In this report, the time-scale separation implicitly assumed by Meyer and Cicolani [4,5] is made explicit 
using forced singular perturbation theory. A new model for the nonlinear controller is derived from these results. 
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The structure of this nonlinear controller is shown in figure 3. Time-scale separation allows the designer to separate 
the aircraft states into groups based on the dynamics of those states. If there is sufficient separation, the designer 
can develop a control law for the slower states while assuming the faster states have no dynamics. That is, if the 
closed-loop plant (inner-loop control system and vehicle) is fast enough to achieve the condition commanded by 
the outer-loop before the outer-loop guidance system computes a new command, the closed-loop plant will be seen 
to have reduced or no dynamics by the outer-loop. This technique provides the designer a simpler system around 
which to build a control system. For example, figure 3 shows the aircraft states separated into very fast states-servo 
actuator states, fast states p, q 9 and r (typically inner-loop states), and slow states V , h 9 a, /?, 9 , and <f> (typically 
outer-loop or guidance states). In fact, there is yet a slower group of states x, and y (used in navigation). The outer- 
loop guidance system must be slow enough to let the plant achieve the commanded states before updating guidance 
commands. The states can be separated into different groups depending on the requirements of the task. 
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Figure 3. Separation of fast and slow time-scale controllers. 

Menon et al [3] provide a detailed derivation of nonlinear flight-test control laws based on the results of Menon 
et al [2] in which the outer-loop flight-test control outputs are fed into an inner-loop CAS. In this latter report, the 
explicit time-scale separation from Menon et al [3] is made implicit, because the CAS is assumed to be much faster 
than the flight-test control. Excellent results were obtained when this flight-test control was tested using a nonlinear 
six-degree-of-freedom simulation containing detailed, full-envelope models of the F-15 aerodynamics, propulsion 
system, and control system. While this work is based on the use of detailed aerodynamic models in the nonlinear 
control, it forms the basis for the application of nonlinear control to the outer-loop control problem addressed by the 
current work. 

A NEW APPROACH TO THE NONLINEAR CONTROL PROBLEM 


It will be demonstrated in the following that the requirement for detailed aerodynamic models in nonlinear control 
laws can be eliminated by using accelerometer and rate gyro outputs. This simplifies the nonlinear controller in 
two ways 
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1. Aerodynamic force calculations are not required in the controller, and 

2. the nonlinear controller to a large extent becomes “model independent” given the same command paths into 
the inner-loop control system. 

Both of these features provide the potential for major design, implementation, and execution time benefits. The 
model-independence feature is perhaps the most significant contribution, because it offers the possibility of designing 
a generic, general-purpose outer-loop control law that is independent of the aircraft to which it will be applied. 
Obviously, complete model independence can only be achieved when one can represent vehicles by identical models. 
For the controller derived and implemented in this paper, only nonlinearities due to aerodynamics and equations of 
motion arc considered. While other nonlinearities exist, they are not accounted for and are a likely cause of some 
of the problems observed in the controller. These are nonlinearities in the inner-loop control system such as stick 
mass inertia, and command scaling. Such nonlinearities contribute to the dependence of the controller on the air- 
craft model. 


Conceptual Structure for Nonlinear Control Law Design 

The nonlinear control laws developed for the FTTC are based on the generic model of the aircraft for outer-loop 
control shown in figure 4 rather than Meyer’s generic aircraft shown in figure 2. The primary difference is that while 
the technique of Meyer et al provides for direct control of surface displacement to control rotation of the aircraft in 
pitch and roll, the technique used in the FTTC relies on the CAS to provide surface displacement controls to achieve 
the commanded values of roll rate and incremental load factor. 



The structure used to develop the nonlinear control laws for the FTTC is shown in figure 5. This structure differs 
from the formal structure for advanced flight-control systems introduced by Meyer [6,7] and shown in figure 1 . The 
structure used for the FTTC uses states based directly on the measured (or reconstructed) aircraft states rather than 
the error between the linear model and the linear system (aircraft and WT-map). The advantages of the current 
approach are that determining the apparent linear system (described later) is fairly straightforward and potential for 
poor tracking due to modeling error is significantly reduced. 
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Command Pseudo- WT-map Aircraft plant 

generator controller 


Figure 5. Structure used for flight-test trajectory control (FTTC). 


A significant difference between Meyer’s formal structure and the one used for the current work is that the 
nonlinear plant in Meyer’s structure is assumed to consist of only nonlinear state equations whereas the current 
work assumes a measurement equation 

z = h(x) ( '0) 

The T-map for the structure shown in figure 5 is defined as 

y = T(z) < n > 


however, most often, T is just the identity matrix. The T-map for Meyer’s controller shown in figure 1 is given by 
equation (4). 


y 

V 


T(x, u) 


The W-map for the nonlinear controller shown in figure 1 is defined by 


u - W(y, v) 


( 12 ) 


Whereas, the W-map for the FTTC is 


u = W( v, z) 


(13) 


While the difference could be reconciled by stating that the measurement function h (x) is included in the T-map, 
there is a subtle difference in that the values in the measurement vector are generally obtained straight from aircraft 
sensors, not computed. 
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Linear and Nonlinear System Formulations 

The plant around which the FTTC is designed includes the aircraft, the CAS and actuators, and the instrumen- 
tation system. This plant is modeled by the equation 

X = fp(x) + gp(x)u 
z = h(x) 

In this model, the number of states must be the same as the number of controls. 

When taken with the WT-map (inverse linearizing transforms) and ignoring the fast dynamics of the CAS and 
other components by applying time scale separation, the system can be rewritten as being linear with the follow- 
ing formulation 

x = Ax + Bv (16) 

This formulation is known as the apparent linear plant. When properly developed, this system is in the Brunovsky 
canonical form. The state vector x includes relevant states selected from the nonlinear system as well as the required 
integral states. 


(14) 

(15) 


Derivation of the Linearizing Transformation 

The input to the apparent linear system, v , is defined as 

x = v (17) 

This equation and equation (14) are then combined so that 

v = fp(x) + gp(x)u (18) 

This is the linearizing transform for state feedback. 

Next, the terms f p (x) and g p (x) in equation (15) must be transformed into terms that are functions of the 
measurement vector, z 


f p (x) = f m (z) (19) 

gp(x)=g m (z) (20) 

So that equation (14) becomes 

X = fm(z) + gm(z)u (21) 

This transformation eliminates the need for aerodynamic models in the linearizing transformations of systems de- 
veloped previously. The linearizing transformation as a function of measurements is given by 

v = f m (z) + g m (z)u (22) 


Each real control input into the plant can then be separately determined by solving equation (22) for each non- 
linear input Uj so that 

Ui = gmi(z) -1 [ Vi - fmi(z)] (23) 

This equation represents the W-map which is the inverse linearizing transformation for the real input to the nonlinear 
system found in figure 5. The T-map in the same figure maps the plant outputs (measurements) to the linear system 
outputs. In this formulation, the T-map is identity. Implicit in this equation is the assumption that g™ ( z) i- 0 . 
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Linear Control Law Design 


In the formulation of the FTTC, the linear control law is a simple tracker as opposed to the explicit model follower 
used by Meyer et al [7]. The linear control law or pseudocontrol receives its input, e , from the feedback, y , of the 
apparent linear system and the reference or command trajectory, y re /» so that 

e = y re /-y 

For the tracker, the first step of the design is to determine the nature of the command to be used (that is, step, 
ramp, or parabola). From this, the designer can determine the number of integrators required for the control law. 
Once this has been established, the pseudocontrol gains can be determined from the apparent linear system using 
modern or classical control techniques. 

Nonlinear System Equations 

For the aircraft described in figure 4, if the inner-loop control states are ignored, the aircraft can be represented 
by a nonlinear system equation in the form of equation (2) 

x(t) = f [x(0] + g[x(t)]u(0 ( 25 ' ) 


with 


x = 


7 

a 

4 > 

e 

h 


(26) 


For a given maneuver, three states to be controlled are selected from these five states. The controls for this system 
are the inputs to the inner-loop controller, and are defined for this application to be 

(27) 



Ul 



u = 

U2 


Pcom 


_ U3 


Vcom 


where 


1 


(28) 


^ ^Zcom 

The nonlinear system state equations used are the standard wind axis equations of motion which correspond to 
the states in the state vector and can be found in reference [34], They are 

1 


V = —[ —Dcos [3 + Y s\n 0 + Xj cos a cos (3 + Zrsinacos/3 

771 

- mg( cos a cos /? sin 9 - sin /? sin <(> cos 9 - sin a cos /? cos <}> cos 9) ] 

a, = 1 f — L + Zt cos a — Xj sin a + mg{ cos a cos (f> cos 6 + sin a sin 9) ] 

Vm cos (3 

+ q — tan (3 ( p cos a + r sin a) 

4> = p + q sin 4> tan 9 + r cos <f> tan 9 
6 = q cos (f> — rsin 4> 

h = — [ sin 9{X T - D cos a + L sin a) - (^) sin <£cos 6 

grn 

-cosmos 9(Zj - Dsin a - Lcosa)] - 1 


(29) 

(30) 

(31) 

(32) 

(33) 
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These equations assume the acceleration of gravity remains constant and side force due to thrust is zero. Of 
these, equations (31) and (32) are already in measurement feedback form. Equations (29), (30), and (33) need 
to be transformed into measurement feedback form. This is done in appendix B, where the linearizing transforms 
are derived. 

If we assume that a maneuver begins at t = t 0 = 0 with the aircraft at the trimmed state, so that x = 0, then 
equation (25) can be written as 

x(0) = 0= f [ x ( 0) ] + g[x(0)]u(0) 

If we define the inputs for the trim condition to be Uj (0) = 0 for i = 1 , 3 , then it follows that 

f[x(0)] = 0 

Extensions to the General Trajectory Control Problem 

Kato and Sugiura [35] present results for the general trajectory problem that, like the results presented herein, are 
based on the work of Meyer et al [4,5,26]. In their article they apply transformations from the wind-reference-frame 
Euler angles to the body-reference-frame Euler angles to show how a general reference trajectory can be controlled 
using nonlinear control. These results are related to those presented by Cicolani and Weissenberger [36]. 

The results of Kato and Sugiura combine the inner- and outer-loop control and transform a general flight trajec- 
tory specified in terms of position variables into vehicle surface commands using the basic nonlinear control methods 
and aerodynamic models. As has been shown in the present paper, the need for these aerodynamic models can be 
eliminated from the outer-loop control law using measurement feedback. It is unclear whether the present results 
can be extended to the inner-loop problem. Further research will address this topic. 

From the point of view of outer-loop control, the problem of navigation requires looking at the aircraft as a 
point in a known inertial system. Although this problem is not addressed directly in this paper, there are no obvious 
reasons why this approach would not work adequately for this purpose. The use of another layer on the outer loop 
may be required, as well as the addition of the inertial position states (x, y , and z ), to the system equations. 

Control Laws 


The control laws for the nonlinear FTTC are driven by error signals obtained by computing the difference be- 
tween the linear state vector and the reference state vector prescribed by the desired aircraft trajectory. The output 
of these linear control laws are pseudocontrols which are transformed, using the W-map, into commands for the 
aircraft-inner-loop system. 

There are three different trajectories contained in the FTTC; a level acceleration maneuver, a pushover— pullup 
maneuver, and an excess-thrust windup turn maneuver. The level acceleration maneuver is a wings level, constant 
altitude maneuver with Mach number increasing (or decreasing) at a specified rate. The pushover-pull up maneuver 
is a wings level, constant throttle maneuver with angle of attack varying at a specified rate from trim to a minimum 
value then to a maximum value then back to trim. The excess-thrust windup turn is a constant altitude, constant 
Mach maneuver with angle of attack increasing at a specified rate to a specified value. 

The control laws used for the implementation discussed later in this paper are described here. A complete 
derivation of each of the control laws shown in the following three subsections is found in appendix B. Basically, 
there are three command loops available to be controlled by the FTTC. These are pitch, which is controlled by 
incremental load factor command, roll, which is commanded by roll rate command, p com ; and thrust, which 

is controlled by throttle command, rj com . This is not to say that these are the only commands for which control laws 
may be developed, but that these are the commands required by the implementation used to validate this approach. 
It would be fairly simple to develop control laws to meet other command requirements. 
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Pitch Control Laws 


The pitch commands for the FTTC are computed using two different control laws and transformation equations. 
The level-acceleration-deceleration maneuver uses an altitude command; the excess-thrust windup turn and the 
pushover-pullup maneuver use an angle-of-attack command. 

The intermediate command of vertical acceleration is computed in the altitude command pitch loop using the 
following linear control law 

h CO m = kh D e h + khpCh + / e hdt ( 34 ) 

•/to 

with 

e h ~ ^Tef ^ 

e/j — h re f h 

Then A is obtained from h com using the inverse linearizing transformation equation (appendix B, eq. (B-28)) 

'hco m ~ ax sin 9 + a, sin (/.cos 6 + 1 (35) 

^ n icom _ + cos Q CQS <f) 

The second pitch control is the angle-of-atlack command loop. The linear control law for this loop is 


Qcom — f 


(36) 


where 


= Qfre/ - « 


From at com an intermediate pitch rate command is computed using the inverse linearizing transformation (appendix B, 
eq. (B-59)) 


qco rn = aicom + tan (3{p cos a + r sin a) 

1 [ — g(a n cos a + a x sin a - cos acos<£cos 6- sin a sin 0)] (37) 

V cos (3 

The pitch rate command q com is then used to determine the A n* com using the following inverse linearizing transfor- 
mation (appendix B, eq. (B-54)) 


A = - 1 + a„ + 


V cos /3 


(9com 9) 


(38) 


g cos a 

These two transformations could be combined into one transformation, however, the pitch rate command generated 
by equation (37) must be filtered because of noise and system delays. This effect is discussed in the Implementation 
Experience section. The separation of these equations is left in here for clarity. 


Roll Control Laws 


The lateral commands for the FTTC are computed using two different control laws and transformations. For the 
level-acceleration-deceleration and pushover-pullup maneuvers, the reference roll attitude <f> Te f is known and 4> com 
is computed using the linear control law 


4* com — 



(39) 
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where 


e 4> ~ tfrref $ 

and p com is computed using the inverse linearizing transformation equation (appendix B, eq. (B-70)) 

Pcom = ficom ~ tan 6( q sin <f> + r cos <f > ) 


( 40 ) 


The control law computation for the excess-thrust windup turn requires the use of two layers of outer-loop 
controllers. Each outer loop must be slower than its inner loop. For this control loop, the outer-loop linear control 
law computes an h com using the same linear control law shown in equation (34). Roll attitude command is then 
computed using the following inverse linearizing transformation equation (appendix B, eq. (B-81)) 


4>ref = ~ COS 


O* sin 9 - h com - 1 


cos 9 




2 + o 2 

y u n 


This intermediate command is used in the linear control law 



(41) 


tficom 

J to 

from which p com is computed using the inverse linearizing transformation equation (40). 


Thrust Control Laws 


Two throttle commands were required for the maneuvers defined. The first simply holds the throttle position 
constant at the current or some predefined position 


Wcom ~ Vref 


The throttle command for all other maneuvers is computed using the linear control law 


where 


Mcom — 


f If 

~ + k Mp e M + / k M] e M dt 

J to 


e/j = M re f - M 


(42) 


eu - M re f — M 

which is used to establish a thrust command T com from the inverse linearizing transformation equation (appendix B 
eq. (B-114)) 


'^' com cos a cos (3 { ma ^ com ~~ a v s ' n P + cos a x cos a — a n sin a) ] — Xj cos a cos /? 

+ mg(cos a cos /?sin 6 - sin ^sin <f> cos 6 - sin a cos 0 cos cos 6) (43) 

The actual throttle setting is determined by using a model of the thrust-throttle relationship (without dynamics) and 
inverting it so that 

Vcom = fj 1 ( h, V, Tcom) (44) 

The thrust term Xj in equation (43) is obtained using the same model with a simple lag to simulate the engine lag. 
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Control Law Coupling 


The inverse transformation described in equation (23) transforms output feedback and linear controls to the 
nonlinear inputs required by the nonlinear plant. This transform provides a decoupled connection between the linear 
control law and the nonlinear plant. However, by decoupling the CAS inputs, errors caused by command limits 
and rate limits can be introduced to the inputs. For instance, if the vertical acceleration guidance law commands a 
load factor greater than the 4.0-p limit, the Mach guidance law linearizing transformation would get an erroneous 
input for normal acceleration. In addition to the limits problem, for decoupling, a different three-axis guidance law 
(with corresponding linearizing transformations) would be required for each maneuver. This adds complexity to the 
programmers task. Layering of outer-loop guidance laws presents another problem because, although the layers for 
each command path are separated in time scale, it is not necessarily the case that there is separation of time scales 
between the inner layer of one path and the next layer of another path. 

The solution chosen for this application is to treat equation (23) as a scalar equation, and to feed back the 
measured values of the terms that are inputs to other guidance paths. This way, limits and other nonlinearities in one 
path do not affect another, although there is some delay associated with this method. This solution also simplifies 
the programmer’s task in that there is only one inverse transform required for each command and maneuvers that 
have identical guidance laws in any one particular channel (such as the angle-of-attack rate guidance law used by 
both the windup turn and pushover-pullup) can use the identical piece of code. 

SYSTEMS DESCRIPTION 


The next three subsections describe the systems and software used to develop, evaluate, and test this approach. 
As has been stated before, the theory and control laws developed in this paper are of generic application. There is 
no requirement to use a specific aircraft or type of aircraft to implement them. The control laws were developed for 
an aircraft which accepts incremental load factor, roll rate, and PLA commands. It should be clear, however, that 
similar control laws could be developed to meet almost any command requirements, because the approach would 
remain the same. The feedback gains selected should, however, keep the aircraft capabilities in mind. Therefore, 
no effort will be made to give an in-depth description of the systems used. Instead, descriptions will be given to 
provide adequate information regarding how the experiment was performed, and to describe drawbacks which had 
or were expected to have noticeable impacts on the results of the experiment. 

Simulation Development System Description 


The system on which the FTTC was developed and evaluated is shown in figure 6. It consists of two computers; 
the simulation computer and the control law computer. The simulation computer contains a six-degree-of-freedom 
nonlinear simulation of the F-15 HIDEC aircraft with full-envelope aerodynamic and thrust models, F-15 CAS 
and actuator models, and atmospheric and sensor models. The simulation computer also contains a simulation of 
the discretization effects and pure delays of the actual flight system. The control law computer contains only the 
FTTC software. The two computers communicate with each other through dual-port shared memory. This system 
is designed to simulate, as accurately as possible, the remotely augmented vehicle (RAV) system. The RAV system, 
on which the controller was flown, is described in the next section. 
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Figure 6. FTTC simulation development and evaluation system. 


Simulation sensor information is taken from the simulation and passed through a simulation of the discretization 
effects and pure delays of the RAV downlink telemetry system before it is passed through shared memory to the 
FTTC. The FTTC uses the simulated sensor information, user-defined maneuver specifications and the control laws 
described earlier to compute commands to be sent back to the simulation. These commands are passed back through 
shared memory and go through a simulation of the uplink telemetry discretization effects and pure delays. The values 
are then passed to the F- 15 CAS model in the simulation computer. The simulation computer integrates the equa- 
tions of motion. Results are then displayed on instruments in a fixed-base cockpit, recorded to disk, or recorded 
on a strip chart. The simulation and controller are monitored and controlled by operator terminals attached to 
each computer. 

This system allows development and evaluation of the FTTC by engineers and pilots. It also provides an effective 
means for software validation, flight system verification and validation, and flight test. 

The simulation is a very good representation of the flight system with two exceptions. First, the simulation runs 
at a rate of 50 Hz. While this is fine for representing the aircraft aerodynamics, it is two and a half times as fast as 
the actual aircraft onboard systems which run at a rate of 20 Hz. However, the simulation in its current configuration 
cannot be operated at 20 Hz. It is possible that aircraft response may vary from simulation response for this reason. 
This was not expected to have a large effect on the flight-test results. 

The second item is sensor noise and bias. While the aircraft sensors are meticulously calibrated, it is impossible 
to eliminate the noise of the sensed data without increasing the loop delay of the overall system. It is also very 
difficult to model the noise and bias of the sensors in the simulation system. This was expected to have a noticeable 
effect on the flight system. 


Flight-Test System Description 

The F-15 HIDEC aircraft and the RAV system shown in figure 7 were used to flight test the control laws. The 
F-15 HIDEC aircraft is a preproduction F-15 with Pratt and Whitney 1128 engines (Pratt and Whitney, West Palm 
Beach, FL), a digital CAS, and a digital electronic engine control (DEEC) computer [1], The RAV system has both 
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downlink and uplink capabilities, and has a ground-based minicomputer in which control laws may be implemented 
for flight test. This system eliminates the need for large computational capabilities onboard the aircraft. 



Figure 7. Remotely augmented vehicle (RAV) system used to flight test the FTTC. 


The RAV system is designed so that software developed and tested on the control law computers of the simulation 
systems may be transported directly and run on identical computers in the RAV system. In fact, it is a requirement 
that no changes be made between the software tested on the simulations and the software used in flight. This enables 
the engineer to perform verification and validation tests cheaply and easily on the simulation before using the aircraft 
and RAV system which have limited availability and cost much more to operate. 

The drawback, however, is that an additional 80 to 100 msec of pure delay are added to the control loop. This de- 
lay added to an aircraft system delay caused by the asynchronous systems of up to three frames (or up to 150 msec), 
can make the control task much more difficult if not impossible. For the FTTC, this was expected to cause a notice- 
able problem since part of the fast controller-the linearizing transform-resides in the ground system. 

The allowable command ranges for the controller were set to be ±150 .0 deg/sec in roll rate, -0.75 to +4.0 g in 
incremental load factor (+0.25 to 5.0 g normal acceleration), and full range for both throttles. 


Flight-Test Trajectory Control Software Program 

The linear control laws and nonlinear inverse transformations described in the Control Laws section were im- 
plemented in FORTRAN in a computer program called FTTC, along with a command generator to produce the 
reference time histories of the states to be tracked. An operator interface to the controller was developed to input 
maneuver specifications and to monitor the FTTC and flight-test system. 

This program is designed to run in the control law computer in the simulation and RAV systems. The basic frame 
time of the program is 12 msec. However, only the maneuver selection logic and the linearizing transformations run 
at this rate. The command generator and pscudocontroller run every other frame. This is done because the command 
generator and the pseudocontroller really represent the outer-loop controller (slow controller in fig. (3)), while the 
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linearizing transformations and the CAS form the inner-loop controller (fast controller in fig. (3)). The linear system 
to be controlled is represented by the linearizing transformations, CAS, and aircraft together. 

The frame time was selected because it was the fastest that the software could run without getting frame overruns. 
A slower frame time would have added to the overall loop delay of the system. However, this means that the 
controller sees constant values of sensed data between downlink data frames which have the same rate as the aircraft 
system (50 msec). This could become a problem during very dynamic maneuvers. 

DISCUSSION 

It is important to note that in the nonlinear control laws defined by equations (34) through (44), none of the equa- 
tions use aerodynamic terms and, hence, none of these control laws require aircraft-specific aerodynamic models. In 
fact, no aircraft-specific models are required except for the simple thrust-throttle mapping shown in equation (44). 
These results are significant in their differences from the results reported by previous authors. In those results, 
expressions for lift, drag, and side force would appear in the equations in place of the accelerometer output terms. 

Thus, the claim of model-independence is shown to be valid under two main assumptions: 

(1) the generic model of the combined aircraft inner-loop control system (fig. (4)) is a reasonable representation of 

the actual system, and 

(2) the time-scale separation of the inner- and outer-loop control laws is a reasonable approximation. 

This somewhat ignores the inverse thrust-throttle mapping. However, it is anticipated that the need for us- 
ing this mapping can be eliminated using an additional feedback loop; this development will be the subject of 
future research. 

During simulation and flight test, controller performance was measured by the ability to keep aircraft parameters 
within specified tolerances. These tolerances arc 

±.01 Mach number, 

±100 ft altitude, and 
±0.3° angle of attack. 

In addition to these, it is required that the controller have smooth transitions from the capture of initial conditions to 
the beginning of the maneuver itself. Also, when parameters to be tracked are ramped, the controller is to maintain 
the tolerances while ramping those parameters smoothly. The large tolerance in angle of attack is due to the noise 
and accuracy problems associated with the measurement of that parameter. 

Simulation Results 

Figure 8 shows the time history from a simulated level-acceleration maneuver from Mach 0.75 to Mach 1.20 at 
25,000 ft. The command rate for this maneuver is M com = 0 .01 sec -1 . The plots show that Mach number increases 
smoothly and that throttle command is smooth throughout the acceleration. Altitude and roll attitude are also held 
relatively constant (less than ±10 ft and ±1.0°, respectively) during the maneuver. This maneuver spans a range of 
dynamic pressure from 309 lb/ft 2 to 793 lb/ft 2 . This represents a large Mach range spanning well more than half 
of the flight envelope. This also represents a range of conditions over which the aircraft can in no way be repre- 
sented linearly. 

Figure 9 shows the time history from a simulated pushover-pullup maneuver at Mach 0.90 and 10,000 ft (825 lb/ft 2 
dynamic pressure). The maneuver is initiated from a wings-level, trim condition. Angle of attack is commanded 
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at a rate of a com - —0.5 deg/sec from trim a first to a = 0° and then at a rate of a CO m - +0 5 deg/sec to a 2.5 . 
During the maneuver, throttle is held constant and roll attitude is commanded to remain wings level. The plots show 
smooth ramping of angle of attack throughout the maneuver. However, angle of attack never reaches the condition 
of 0° because the commanded incremental load factor reaches its lower limit before that condition is attained. On 
the other hand, the angle of attack never reaches the desired maximum condition of 2.5° because the integrator in 
the angle-of-attack pseudocontrol has wound up while the commanded incremental load factor was limited. By the 
time the command generator starts the pushover from the maximum angle-of-attack condition, the aircraft has not 
yet achieved 2°. Neither the proportional nor the integral part of the linear control law can act fast enough or have 
great enough magnitude to eliminate this error. The reason is that the CAS does not provide an exact linear control of 
measured load factor based on commanded incremental load factor. This is discussed further in the Implementation 

Experience section. 

Figure 10 shows the time history of a simulated excess-thrust windup turn at Mach 0.65 and 25,000 ft, which 
corresponds to a dynamic pressure of 232 lb/ft 2 . This maneuver begins at the straight-and-level-trim condition and 
increases angle of attack to 12° at a rate of 1 deg/sec. Altitude is maintained to within ±60 ft. Mach number is held 
to within ±.02 until the end of the rollout back to trim. This error, which is in excess of the specified tolerance, is 
caused by a 1.5-sec delay built into the throttle command to eliminate the problem of the throttle toggling between 
core power and afterburner. It is expected that if this delay were not required, Mach number would have been 
maintained within the tolerance. For this maneuver, the aircraft longitudinal, lateral, and thrust controls are coupled 
because of the relatively high angle of attack and steep roll attitude. But while the controller handles this level 
of coupling adequately, there are indications that problems may exist at slower speeds or higher angles of attack. 
Each time the angle-of-attack control changes between a steady-state control and a ramp control, the roll controller 
has difficulty controlling altitude. However, up to this level of coupling, the controller’s performance during this 
maneuver demonstrates that this approach handles this problem. 

Figure 11 shows the time history from a simulated excess-thrust windup turn at Mach 1.2 and 25,000 ft. This 
maneuver reaches an angle of attack of 5°. The plots show that altitude is held to within 40 ft and Mach to within 
Mach 0.005. For this maneuver, dynamic pressure is 793 lb/ft 2 . The performance of the controller during these two 
windup turns proves that the controller is capable of fairly precise control over a large range of dynamic pressures 
(232 lb/ft 2 to 793 lb/ft 2 ) without changing control law gains. 

In fact, simulation evaluation was conducted across a Mach number range from 0.60 to 1 .20 and across an altitude 
range from 10,000 ft to 40,000 ft, with acceptable results except at the high and slow conditions. While these results 
did not fall within the specified tolerances, the controller was still able to perform the maneuver. It is expected 
that the slower, higher flight conditions which have still lower dynamic pressures may be a problem because the 
aircraft dynamics get slow enough to invalidate the assumption of time-scale separation between the inner- and outer- 

loop controls. 

Simulation tests were also conducted with varying amounts of pure delay in the uplink and downlink simulation. 
Total system delays ranging from 120 to 240 msec were examined over a matrix of maneuvers simulated at Mach 0.75 
and 1.20 and at altitudes of 10,000 ft, 20,000 ft, 30,000 ft, and 40,000 ft. The controller performed acceptably at all 
conditions, with the same exceptions indicated in the previous paragraph. 
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Flight- Test Results 


The FTTC was flown between June and August, 1989. During these flights, pseudocontroller gains were adjusted 
and minor changes were made in the control law to improve performance of the controller. These gain changes 
were required primarily because of nonlinearities in the aircraft CAS that were not accounted for in the controller 
and because of measurement noise. The controller did fulfill the expectations of providing adequate control of the 
aircraft over a large range of flight conditions, although some work was still required to get the quality of results 
expected. Four of the maneuvers flown are presented in figures 12 through 15 and are discussed here. 

The first maneuver, which is shown in figure 12, is a level acceleration maneuver from Mach 0.70 to 1.20 at a 
rate of 0.01 Mach number per second. The roll attitude is held relatively constant with an error of less than 2.0°. 
Mach number increases at a constant rate except for a short time around Mach 1.0 (at approximately the 30- sec 
mark). This effect is caused by sensor error as the shock wave forms at the tip of the nose boom. 

A much larger transonic effect can be seen in the altitude, altitude rate, and commanded incremental load factor 
plots. The cause of phenomena seen in these plots can be traced to the source of these measurements. The altitude 
and altitude rate measurements are taken from the inertial navigation system on the F-15 aircraft. Inertial altitude 
rate is derived by filtering and integrating a set of accelerometer and rate gyro inputs. Therefore, it is not affected by 
the shock wave across the pitot-static system on the airplane. However, inertial altitude is obtained by integrating 
the value of altitude rate and removing accumulated biases by comparing inertial altitude with the altitude from the 
static pressure system while the aircraft is in nearly unaccelerated flight (that is, level flight). Hence, as the aircraft 
passes through the transonic regime, the altitude tends to decrease because of the sensed increase in static pressure 
before the shock. The inertial system filters this information into its estimate of altitude. The resulting value is sent 
down to the controller which attempts to correct by climbing. The reverse occurs as the aircraft achieves supersonic 
flight and the static system is behind the shock wave. This is the cause of the apparent discrepancy between altitude 
and altitude rate in the time histories. However, the inertial system is slow in filtering the pressure altitude into the 
inertial altitude, so the effect, while noticeable, is fairly small and leads to an altitude error of approximately ±60 ft. 

The pushover-pull up maneuver shown in figure 13 was flown at Mach 0.90 and 10,000 ft with angle of attack 
varying between 0° and 2.5° at a rate of 0.5 deg/sec. Angle of attack is tracked well in this maneuver although sensor 
noise does affect the commanded value noticeably. From these time histories, it is difficult to determine whether the 
noise in the sensed angle-of-attack trace (which is unfiltered) is caused by noise in the commanded incremental load 
factor trace, or vice versa. Other maneuvers were performed in which commanded incremental load factor remains 
constant yet there is still noise on the angle-of-attack trace. The problem of noise on angle of attack shows up in the 
next two maneuvers and will be discussed later in more depth. 

Figures 14 and 15 show excess-thrust windup turns. The first one is flown at Mach 0.65 and 25,000 ft, while the 
second one is flown at Mach 1 .20 and 25,000 ft. The maneuver at Mach 0.65 tracks altitude poorly with a maximum 
error of approximately 150 ft before it starts to roll out to straight and level. This maneuver also has difficulty 
tracking Mach number with an error of approximately ± 02 Mach for the same reason as the simulated maneuver 
discussed in the last section. The level of control coupling caused by the angle of attack and steep roll attitude 
presents a greater problem to the controller in the flight system than in the simulation system. This coupling along 
with the throttle delay contributes to the large altitude excursion. 

The second excess-thrust windup turn, which is performed at Mach 1.20, does not show the pitch and thrust 
coupling problem. In fact, altitude is kept to within ±25 ft, and Mach number is kept to within ±.008 of the 
commanded value. Here, however, angle of attack is not ramped smoothly to the desired value of 5°, and even 
though the commanded incremental load factor is at its maximum value, the angle of attack measured is just over 3°. 
This problem and a similar one between commanded and measured roll rate will be discussed in the Implementation 
Experience section. 
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The roll racheting effect was the effect during flight test which the pilots found least acceptable, and in some cases 
unacceptable. This effect can be seen at the beginning of the roll into the turn. Roll racheting varied in magnitude 
and frequency depending on the conditions of the turn. In the slower or higher turns the effect was not noticeable to 
the pilot. The origin of this effect is unclear. While it did not show up in the simulation evaluation of the system as 
flown, it could be made to appear in several ways by varying gains, maneuver requirements, delays, etc. Therefore, 
no definite source was determined. This effect is also discussed further in the Implementation Experience section. 

The origin of the noise in the angle of attack and commanded incremental load factor seen in the pushover-pullup 
maneuver can be explained, at least in part, by the maneuvers in figures 14 and 15. During both maneuvers, the com- 
manded incremental load factor reaches its maximum limit and is held constant there for a short time. However, noise 
in the angle-of-attack trace is significant (up to approximately 0.3°) and is probably a result of either measurement 
noise, turbulence, or both. Angle of attack provides the feedback for the intermediate command of angle-of-attack 
rate (see eq. (36)), as well as heavily influencing the inverse transformation shown in equations (37) and (38). The 
commanded incremental load factor computation seems to amplify the noise on angle of attack. While this is not 
a divergent process, filtering sensed angle of attack will be necessary to provide smoother tracking of commanded 
angle of attack. 

The results of the FTTC flight test indicate that the approach using linearizing transforms with measurement feed- 
back is a valid approach for developing outer-loop controllers for nonlinear systems. Tracking reference trajectories 
is acceptable over a large range of altitudes and Mach numbers using only a single set of gains. 

In general, the flight-test data matched the simulation data with the exception of sensor noise, which was not 
modeled in the simulation. The problem of roll ratcheting warrants further research. 
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Figure 13. Flight data of a pushover-pull up maneuver flown at 10,000 ft, Mach 0.90 with angle of attack varying 
between a = 0° and a = 2.5°. 
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Implementation Experience 


There were several factors in the implementation that caused problems to the experiment. These factors are 
discussed here to provide some insight to problems which may be encountered in future implementations. 

During the development, simulation, and flight test of the FTTC, some adjustments in gain values were required, 
although the objective was to develop a technique which would eliminate gain adjustments. At least part of the reason 
for this is the inaccuracy of the modeling of nonlinearities in the inverse transformations. While the transformations 
used account for the nonlinearities in the equations of motion and in the aerodynamics of the aircraft, they do not 
account for the nonlinearities in the aircraft control system. One example of this is discussed later in this section when 
the roll ratcheting problem is examined. Further research is required to determine if modeling these nonlinearities 
actually minimizes gain adjustments. 

* n lhis implementation, the F- 15 CAS did not control the aircraft to provide the roll rate commanded by the 
FTTC. In fact, roll rate command was two and a half times as great as the achieved roll rate, as can be seen by the 
flight data from the turning maneuvers. Also, flight data indicate that this scaling factor is not constant, but varies 
slightly with flight condition. This might provide a partial explanation for the roll ratcheting problem discussed in the 
Flight-Test Results section, since it could cause the loop gain to be too high at some flight conditions. Additionally, 
the sensed and commanded roll rates from flight-test data do not compare to those from the simulation in terms 
of their relative magnitude. This points out that there is a modeling inaccuracy in the simulation. Such modeling 
inaccuracies explain the need to adjust gains once flight test is underway. 

Flight data also indicate a similar scaling factor problem between incremental load factor and achieved load 
factor. This presents a problem in accurately tracking the reference parameter in the pitch axis, although there 
is no similar ratcheting problem in pitch. This along with the limited range in incremental load factor command 
significantly limits the capability of the controller to perform maneuvers. 

As discussed in the Control Laws section, the same linear control law (eq. (34)) was used for altitude tracking 
in both the pitch and roll control loops. However, the actual implementation of these two control laws required 
different gains for each of the loops. While it is not obvious that this should be the case, examination of the control 
loops provides insight as to why this is necessary. The h com control law of the altitude commanded pitch control loop 
is airly fast since it has only one outer loop. The h com control law of the altitude commanded roll control loop is 
designed to be slower since it is the outermost loop in a multilayered control loop and must observe the requirement 
of time-scale separation keeping it slower than any of the inner loops of this control. 

The implementation of the control laws in the FTTC does not provide control decoupling for the reasons de- 
scribed in the Control Law Coupling section. As can be seen in the simulation and flight data of the maneuvers, 
this does not appear to be a problem. Very good performance is achieved in highly coupled maneuvers despite the 
decoupled inverse transformation equations. Coding is much simpler. Verification and validation testing is much 
less complicated and provides a higher level of confidence in the software. 

As was discussed in the Flight- Test Results and Pitch Control Law sections earlier, one of the contributors to 
t e noise in the angle-of-attack based incremental load factor command is the sensitivity of the intermediate q com 
command to noise in measured angle of attack. The solution used here was to put q com through a low-pass filter. 

e difference between q com and measured q was then computed, and run through a proportional-integral control 
aw to eliminate biases and standoff errors. The output of the proportional-integral control law would then be used 
in equation (38). This solution helped a great deal in eliminating high-frequency content of the command, and in 
improving tracking of the reference angle of attack. Another solution that was not tried would be to filter measured 
angle of attack. While this would eliminate the high-frequency content of the angle-of-attack signal, it would also 
add delay to the control loop. This was not tried because of the added delay. 

One interesting unforeseen problem uncovered in the attempt to find the source of roll ratcheting was uncom- 
manded stick movement. When a large step was input into the roll channel of the CAS in a turn, the aircraft rolled 
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sharply into the turn. The stick, because of its own inertia, would tend to resist the rolling motion, and put a roll 
command into the CAS in the opposite direction. This could cause the roll ratcheting effect experienced by the pilots 
and seen in the data. Maneuver conditions requiring less sharp changes in roll rate were attempted which reduced 
the ratcheting slightly. This effect and the command scaling effect in the roll and pitch command paths should be 
accounted for in the inverse transformation. However, it would be preferable to eliminate feedback through the 
control stick by setting stick commands to zero in the CAS while the system is engaged. 

Comparison of Simulation and Flight-Test Results 

The conditions of the simulation maneuvers and of the flight-test maneuvers (figs. (8) to (15)) were selected to 
provide a direct comparison between simulation and flight-test data. The only discrepancy is that the simulated level 
acceleration begins at Mach 0.70 while the level acceleration performed in flight starts at Mach 0.75. 

While in general the simulation and flight-test results agree, it would clearly enhance the ability to analyze 
performance and select gains if the simulation were to run at the same rate that the aircraft systems run, and if the 
FTTC were to run at the same rate as the simulation or slower. In addition, modeling airdata noise in the simulation 
would have enabled the designer to fine tune the gains and the control laws for the flight system. 

The level-acceleration maneuvers compare very well except for the transition through the sonic region in the 
flight data, which was discussed in the Flight-Test Results section. The transient at the end of the maneuver, as the 
controller flies the aircraft from the final condition of the maneuver to straight and level, unaccelerated flight, is 
much sharper in the flight data than during the simulation maneuver. The cause of the transient is the reset done in 
the pseudocontrollers when the controller switches from the maneuvering phase to the exit phase. During the exit 
phase, the controller recovers the aircraft to straight and level, unaccelerated flight. This reset is done to eliminate 
accumulated biases in the pseudocontrol laws enabling the controller to drive the aircraft to straight and level flight 
without first unwinding an integrator. The reason for the difference in transient magnitude between the simulation 
and flight data is that the aircraft has not recovered completely from the erroneous altitude data received while flying 
through the transonic region. 

One other discrepancy can be seen in the incremental load factor command for each maneuver. For the simulated 
maneuver, the controller maintains zero or slightly negative values, while in flight test, the average value tends to be 
about 0.3 g. The probable cause for this is that the pilot had trimmed the aircraft slightly nose down and the controller 
was trying to eliminate this bias (which would be taken out by the integrator in the altitude pseudocontroller, causing 
the reset transient discussed in the previous paragraph). 

The simulated windup turn at Mach 0.65 compares well with the flight tested turn. Mach and commanded PLA 
are nearly identical. However, there is no high-frequency content in the commanded roll rate for the simulated 
maneuver as there is in flight. This is probably because of noise in the measurements. In addition, while both 
trajectories achieve the target angle of attack, the simulated trajectory does not require nearly as much commanded 
incremental load factor. There is no clear explanation for this at this time. 

Finally, the simulated and flight-tested windup turns at Mach 1 .20 match rather well. The main difference is that 
the roll ratcheting seen in the flown maneuver is not seen in the simulated maneuver. Even so, the roll attitudes for 
the two maneuvers are similar, and the altitude time histories are also close. The differences in Mach number and 
commanded PLA are small and can be attributed to the fact that while the test-flown maneuver starts slightly above 
the specified initial Mach number, the simulated maneuver starts a little below the specified initial condition. The 
angle of attack and commanded incremental load factor are similar when noise and a slight difference in initial trim 
condition are considered. 


33 


CONCLUDING REMARKS 


Research has been ongoing at the Dryden Right Research Facility to define a technique for designing flight- 
test trajectory controllers. Many techniques have been investigated. The most promising technique of the previous 
research used linearizing transformations using state feedback. 

In this paper we have presented a new approach to the outer-loop trajectory control problem for developing non- 
linear control laws using measurement feedback. A controller was designed which applies this approach specifically 
to the problem of flight-test trajectory control. Testing was done on a full-envelope, nonlinear, si x-degree-of- freedom 
simulation and in flight using an F-15 aircraft. Right-test results are presented. 

Simulation and flight-test results show that the flight-test trajectory controller can control an F-15 aircraft through 
level accelerations, pushover— pullups and windup turns over a large portion of the flight envelope. Some problems 
are encountered at high and slow (low dynamic pressure) conditions where time-scale separation between the fast 
and slow controllers disappears as aircraft performance decreases. 

The approach combining inverse linearizing transformations with measurement feedback for designing outer- 
loop controllers for high performance aircraft has been validated by the results in simulation and flight testing. Since 
this approach does not require implicit or explicit aircraft models in the design (with the exception of a simplified 
engine model), it may be easily applied to other aircraft requiring similar inputs to the control system. The limitations 
of this approach fall into three categories 

• Separation of time scales must be maintained, limiting the overall performance of the combined control and 
guidance system, 

• nonlinearities found in the apparent linear plant must be accurately modeled so that the controller can be 
designed to deal with them over the entire envelope, and 

• controller gains selected during analysis will probably require adjustment because of modeling inaccuracies 
in the analysis system, or because of nonlinearities which are not accounted for in the design of the inverse 
linearizing transformations. 

The first limitation implies that an integrated approach can provide higher dynamic response because no time- 
scale separation is required. Since the time scales must be separated sufficiently for worst case plant dynamics, 
a penalty will be incurred when plant dynamics get faster. This may justify the use of gain scheduling for the 
linear controllers. 

The second limitation is illustrated by the roll acceleration fed back through the stick which is described in 
the Implementation Experience section of this paper. While the inverse transformations as derived in this paper 
account for all nonlinearities in the equations of motion and the aerodynamic characteristics of the airplane, there 
are nonlinearities in the CAS which remain unaccounted for, and cause problems for the controller. Also, the effects 
of the engine nonlinearities are not completely understood at this time, however, they seem to be relatively minor. 

The third limitation is common to all outer-loop design techniques. It is not unusual for control systems to be 
modified because of response in the simulation or in flight. 

Dryden Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, June 1, 1990 
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APPENDIX A 

ILLUSTRATIVE EXAMPLE OF NONLINEAR CONTROL 

To illustrate the nonlinear control method of Meyer et al, we present an example in which we consider a two-state 
model of an aircraft flying a level acceleration-deceleration 

(A-l) 
(A-2) 

In this example, V and h are the state variables; sin 7 and tj are control variables. We also assume a simple thrust 
and drag model defined about a trim condition represented by the subscript 0 , with 

T = T(V,h)-T 0 ( A ' 3 > 


•, rfT - D 

V = 9 sin 7 

m 

h = V sin 7 


and 

D = D o(V -Vo ) 2 

For level acceleration-deceleration, we define the trajectory 


Vref = Vo + Vr'fA t 


( A-4) 
(A-5) 


and 

h re f = ho (A- 6 ) 

The nonlinear system represented by equations (A-l) and (A-2) has the form required by equation (2), with 


X = 

X X 
■ 

u = 

I 1 

c c 

to — 


' _D 

f = 

m 

0 

1 

• T ' 

gl = 

m 

0 _ 

g2 = 

-9 

V 


V 

h 


V 

sin 7 


(A-7) 
(A- 8 ) 
(A-9) 
(A- 10) 
(A- 11) 


which satisfies f[x(0)] = 0 if x(0) = x 0 , the trim condition. While this is a valid representation of the nonlinear 
system in the form required by equation ( 1 ), it should be noted that this representation is not in the canonical form 
described by Su et al [32], 


We consider the tracking problem by letting 


yi 


1 

1 

l 

. y2 . 


, 

1 

v> 

fc* 

1 


( A-12) 


We will transform this system to the form shown in equation (5) by letting 


y(t) = v(i) 
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so that 


V = 

Vl 


1 

n 

■ 1 


Vref + 9 sin 7 


v 2 


-h 

- 


-V sin 7 


Assuming a linear-quadratic design, this results in pseudocontrols defined by 


Vl 


~ Kn 

Kn ' 

yi 


' Kn 

K n ' 

' Vref - V ' 

V 2 


K21 

K22 _ 

. y2 j 


K 2 \ 

K22 

Kef ~ h 


From equations (A- 13) and (A- 14), the real controls u can be computed so that 


or 


and 


sin <7 = 


7 = sin' 


-folCKe/- V) - Knihref-h) 

V 

zKldVref ~ V) - KlAhref - ft) 
7 


’l ~ "yT [^11 (V re f - V) + KniKef - h) - g sin 7 - 7 re /] + — 
which also can be expressed as 

( Vre/ ~ V) [ VKu + 9K2I ] + (hr ef ~ /l) [ VK \2 + 9 K2I ] 


-TO 


V = 


7 


— Vre/ 


+ 

r 


(A-13) 

(A-14) 

( A- 1 5) 
(A-16) 
(A-17) 

(A-18) 


Three features can be observed in this simple example. The first, which is somewhat obscured by the simplifying 
assumptions, is that a single control law is designed to control state variables based on a reference trajectory, that is, 
the inner- and outer-loop control functions are combined. The second and most important feature, apparent even in 
this simple example is that the actual controls depend on the use of the aerodynamic and thrust models represented 
by D and T. Third, although the altitude and velocity command paths are clearly coupled, the control law provides 
decoupled control of the aircraft flightpath angle and throttle position. 
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APPENDIX B 

DERIVATION OF NONLINEAR CONTROL LAWS 
FOR FLIGHT-TEST TRAJECTORY CONTROL 

Nonlinear System Equations 


The primary equations used to derive the WT-maps for transforming the nonlinear systems into the linear repre- 
sentations are given by the following equations from Duke et al [34] 


h = a Xi k sin 8 - a^k sin <j> cos 8 - a z ,k cos <j> cos 8 
w cos a — u sin a 
** ~ V cos p 

<j) - p + q sin 4> tan 8 + r cos <f> tan 8 
V = it cos a cos /? + v sin (3 + ui sin a cos /? 

= _L[ -Dcos P + Y sin f3+ Xt cos a cos f3+ Zt sin a; cos (3 

m 

— mg( cos a cos (3 sin 8 — sin sin <f> cos 8 — sin a cos /3 cos pcos 0)] 

Additional equations used in deriving this mapping are provided by 

h = V ( cos a cos 13 sin 8 - sin p sin <j> cos 8 - sin or cos p cos <f> cos 8) 

a x k = — (Xt - D cos a + L sin a — gm sin 8) 

' gm 

a v k = — ( Y + gm sin $ cos 8) 
gm 

Qzk = — (Z T — D sin a — L cos a + gmcos<pcosd) 
gm 

= Qx sin 0 
O y k = a y + s * n <£ cos ® 
k + cos 0cos 6 

d n — CL Z 

M- V - 

a 


u = 


_ X T - gm sin fl - D cos a + L sin a + ^ sin p _ qV s j n ^cos /? 
m 


V _ gmsin<ftcosfl+_Y + sjn aC0S ^ _ r y C os <*cos p 

m 

Zt + gm cos <£ cos 0 — £> sin a; — L cos a 


tu = 


m 


(B-l) 

(B-2) 

(B-3) 


(B-4) 

(B-5) 

(B-6) 

(B-7) 

(B-8) 
(B-9) 
(B-10) 
(B-l 1 ) 
(B-12) 

(B-13) 
(B-l 4) 
(B-l 5) 
(B-l 6) 


+ q V cos a cos (3 - pV sin /3 

Derivation of Pitch Control Inverse Transformation Equations 

The equations related to longitudinal control, that is, the equations that ultimately result in a A command to 
the generic high-performance aircraft model (fig. (2)), have different forms depending on the maneuver for which the 
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controller is intended. In this section, the derivation of the WT-map for each maneuver is presented. The derivation 
of the apparent linear system equations is shown for the altitude commanded control loop. The apparent linear 
system for each of the other command variables and control variables is similar and will not be shown explicitly. 


Altitude Command 


For the altitude command control loop, we use the apparent linear system defined by equation (16) 

x = Ax + B\ 


with 


such that 


and 


y = Cx 



. V1 1 


’ h ' 


' 1 

0 ' 


' h ' 

h 


0 

1 


h 


h 

h 


0 1 
0 0 


h 

h 


0 

1 


( B-17) 
(B-18) 
(B-19) 

(B-20) 

( B-21) 


The W-map from v into u is defined by 

u = W(v,z) ( B-22) 

For this control law, the linear system input, v, is selected to be 

v = [vi] = [h COTO ] (B-23) 

and the nonlinear system input, u, for the F-15 CAS is given by 

U = An Zco m = Orw - 1 (B-24) 

To derive the linearizing transformation, the measurement equation for h shown in equation (B-l) as a function of 
kinematic accelerations is converted to an equation for h as a function of accelerometer outputs using equations (B-9) 
through (B-l 2) to get 

h = a x sin 6 - a y sin <f> cos 6 + a n cos <J> cos 9 - 1 ( B-25) 

Converting this equation from a measurement equation to a command equation with o„ = a^ m and h = h com gives 

hcom = a x sin 6 - a y sin <j> cos 9 + cos </> cos 6 - 1 ( B-26) 

The linearizing transformation is obtained from this equation and equations (B-23) and (B-24) to get 

iM = a x sin 6 - a y sin <f> cos 0 + ( u + 1) cos cos 6 - 1 ( B-27) 
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The inverse linearizing transformation is obtained by solving the equation above for the real nonlinear input, u, 

so 11131 _ vi - Ox sin & + ay sin eft cos 0 + 1 _ ^ (B-28) 

u ~ cos 6 cos <f> 

The tracker, or pseudocontrollcr, for this control law will track reference altitude rate, h ref , and reference altitude, 
h ef . Reference altitude rate is ramped in and out to achieve changes in altitude in the FTTC application. e 
resultant reference altitude input will then be parabolic which requires two integrators. A good compromise here is 
to use a proportional-integral-derivative (PID) controller around altitude. This will give some damping with the h 
term, as well as some ability to eliminate steady-state bias with the f h term, while not slowing down the control 
too much by having a double integral on altitude. 

For this case, then, the tracking error is computed by equation (24) 


c = yr ef - y 


with 


e = 


eh 

e h 


y ref 


href 

href 


y = 


h 

h 


so that 


eh 


Kef 


' h ' 

. . 


V 

1 


h 


The term eh is integrated once to get / eh which is combined into e to form e so that 


e = 


I eh 
eh 

ei 


The pseudocontrol can then be computed from 

v = Ke = [ kh, kh P kh D ] 


J eh 
eh 

ei 


(B-29) 

(B-30) 

(B-31) 


(B-32) 


(B-33) 


(B-34) 


which can also be expressed as 

vi = kh D e h + kh P eh + kh, f ehdt 

Jto 

With these two formulations, the gain matrix K can be determined from any appropriate modern or classical 
control technique. Again, the designer should be careful when selecting gains so as not to make the dynamics o t e 
pseudocontroller so fast that it can affect or be affected by the dynamics of the inner-loop control. 
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Normal Acceleration Command 


A linear control law and inverse transformation equation are derived here for a pitch command based on normal 

acceleration. While it was not used for any maneuvers flown in flight test, it was implemented and tested in the 
simulation with good results. 

For this command loop, a„ is measured directly on the aircraft so that in the nonlinear system the feedback v 
is simply ’ J ' 

y = a n (B-36) 

It is desirable to track the reference normal acceleration directly with 


A ti ? = a 1 

£ com 'Ve/ 

and 


a ncom ” a n re f ~ O n 


SO 


A = a>w. + o„ - 1 

Then if the pseudocontrol, v 2 , is defined as 

V 2 = Or w 

and the real nonlinear control is defined as it is in equation (B-24) 

U = An *c om 

The form of the linearizing transform equation would then become 


v 2 = u — a n + 1 


the inverse linearizing transformation is given by 

u = i >2 + a n — 1 

The tracking error is defined by 

e a» = a n rcf — Cln 

The error is integrated once and the combined error vector is 


e = 


/ e a„ 


(B-37) 
(B-38) 
( B-39) 


(B-40) 


( B-41) 


( B-42) 


(B-43) 


Since the plant dynamics are algebraic, a proportional plus integral control law is used for this control task. The 
pseudocontrol, then, is computed using the following equation 


v = Kk = [ fc OB/ k anp ] 


I e a n 


( B-44) 
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Angle-of-Attack Command 

The angle-of-attack command control law requires two outer-loop layers as opposed to the single layer outer- 
loop laws required for altitude control and normal acceleration control. The reason for this will be apparent as the 
linearizing transformation for the outermost layer is derived. 

For the outermost loop of this control law, the pseudocontrol, v 3 , is defined to be 

— a ( B-45) 

^3 _ u com 

and the real nonlinear input is shown in equation (B-24) 

« = An z«, m = ~ 1 


The linearizing transform equation is based on the observation equation for angle-of-attack rate given by equa- 
tion (B-2) 

w cos a — u sin a 
a= V cos /3 

The terms w and u in this equation are transformed from the force equations shown in equations (B-16) and (B- 14) 
to kinematic acceleration equations by using equations (B-8) and (B-6) to get 

w = ga Zi k + V(q cos acos - psin /3) (B-46) 

u = ffai.k + y(rsin (3 - qf sin orcos/?) (B-47) 

These equations are transformed into functions of accelerometers by using equations (B-ll), (B-12), and (B-9) 
to get 

w = -g(a n — cos cos 6) + V( q cos a cos (3 - p sin (3) (B-48) 

u = g(a x - sin 9) + 7(rsin p - q sin acosy3) (B-49) 

Substituting w using equation (B-48) and ti using equation (B-49) into equation (B-2) gives an expression for d in 
terms of accelerometer outputs 


[ — g(a n cos a + o x sin a — cos acos <£cos 6 — sin asin 9) ] 

V cos (3 

+ q — tan f3( p cos a + r sin a) 


(B-50) 


While this equation does describe d as a function of o„, taking the derivative gives 


da _ — g cos a (B-51) 

da n ~ V cos (3 

Unfortunately, this derivative shows the wrong sign indicating that as d increases, a n decreases. This would produce 
aircraft motion opposite to what is desired when commanding angle-of-attack rate. 

It is required then to break up this controller into two control laws. The innermost layer of the outer-loop 
controller will make a pitch rate command available to the outermost layer. The outermost layer then will use a 
command equation similar to the d equation shown in equation (B-50). 

For the innermost layer of the outer-loop controller, the linear input is 
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and the nonlinear system input is the same as in equation (B-24) 


u = A ^ mm = “tw - 1 

The pitch rate command equation is obtained by solving the a equation (eq. (B-50)) for a n . The equation is then 
changed into a command equation by replacing a n with and q with q c^. This gives a new equation 


Hom = — — {V r cos/?[g com -d-tan0(pcosa+ rsin a)] 
g cos a 


+ g( cos a cos <j> cos 9 + sin a sin 6 — a x g sin a) } 


Then substituting in equation (B-50) for a gets 


V cos f) 

~ an+( gCOS a HqC<m ~ q) 


(B-53) 


(B-54) 


Finally, after substituting in the definitions of the linear input and the nonlinear real input, V 3 , and u respectively, the 
inverse linearizing transformation becomes 


Vcos/3 

U 1 + O n + ( ) U3 

g cos a 


(B-55) 


Since the input to the a controller is expected to be a ramp, the input into this control should also be a ramp. 
Therefore, only one integrator will be required, and the vectors e and e will be 


® — [fig] — qcom q 


(B-56) 


e = f Cq 


( B-57) 


The pseudocontrol, vj t , can then be found from 


•J, = [ K V K V ] Je 

J C/T 


(B-58) 


The linearizing transform of the outermost layer of this control is found by first solving equation (B-50) for q 
and then replacing the q and a terms with q com and a com . This gives the equation 


Qcom = oi com + tan /?( p cos a + r sin a) 

— — - [ -g(a„ cos a + a x sin a — cos a cos <£cos 0 — sin a sin 0)] 

v cos p 


(B-59) 


Then by defining the linear system input and the real nonlinear input as 


u 3 = a c 


(B-60) 


U3, = Qcom 

and applying these substitutions to equation (B-59) inverse linearizing transform is obtained as 
«3, = i >3 + tan /?(pc os a + rsin a) 

~ V cos ~3 ^ ~y( an cos a + 0l s ’ n a ~ cos aC0S <^ cos 0 — sin a sin 0)] 


(B-61) 


(B-62) 
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It may seem disconcerting that equation (B-50) is used so often. In fact, that equation, using measurement values, 
is used only once in obtaining equation (B-54). In all other instances, equation (B-50) is first transformed into 
a command equation by replacing a with a com , or o„ with a or q with q com . In these circumstances, the a 
equation is no longer a measurement equation. 

The pseudocontrol law for the outer layer will contain only one integrator since the angle of attack is expected 
to be a ramped command. Therefore, the e and e vectors are defined to be 

e = [e a ] = otref ~ (B-63) 


and 



The pseudocontrol, t> 3 , can then be found from 



(B-64) 


(B-65) 


The designer must be aware that gains need to be selected carefully for this control loop since die dynamics of 
the outermost layer of this control law must be slower than the innermost layer. The dynamics of the innermost layer, 
in turn, should be slower than the plant dynamics. Violating either of these requirements will lead to poor control of 
the system or even drive the system unstable. 

After completing this derivadon, the designer should notice that this control law provides not only a method of 
tracking an angle-of-attack command, but also a method of tracking a pitch rate command. 


Derivation of Roll Control Inverse Transformation Equations 


The equations used to provide roll control are derived here. As with the pitch commands, there are several modes 
depending on the reference parameter to be tracked. The real nonlinear command required by the F-15 CAS is roll 

rate, Pcom- That is 


Roll Attitude Command 

For the roll attitude command control law, the pseudocontrol, v 4 , is defined as 

V 4 = 4>com 

The linearizing transformation equation is derived by using equation (B-3) 

</> = p + q sin 4> tan 6 + r cos <t> tan 9 

and substituting in (p com = tj> and p com = p to get 

4>com = Pcom + q sin <f> tan e + r cos <f> tan d 
Incorporating equations (B- 66 ) and (B-67) then gives the linearizing transform 

V 4 = u + q sin 4> tan 9 + r cos <f> tan 6 


(B-67) 


(B- 68 ) 


(B-69) 
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(B-70) 


From this, the inverse linearizing transform is obtained which is 

u = V 4 - q sin <f> tan 0 - r cos <£tan 6 


The pseudocontrol is computed using a proportional-integral control law since most inputs are expected to be 
ramps or steps. For this control law, the e and e vectors are 


e = i e 4>] = 4 >t ef~4> 


and 


e = 


€ 4 > 


Therefore, the equation to compute the pseudocontrol is 

K^j K^ p j 


V 4 


5 H 

H 


(B-71) 

(B-72) 


(B-73) 


The same requirements apply here as in the pitch control laws regarding the speed of the dynamics in the plant and 
in the outer-loop controller. The only difference is that the aircraft is very fast in the roll axis, so that the linear 
controller dynamics can be very fast as well. 


Altitude Command 


The altitude commanded roll control law is another multilayer control like the angle-of-attack commanded pitch 
control. The outer layer of this control law generates a roll attitude reference based on the altitude and altitude rate 
commands. The reference roll attitude is then used in a roll attitude command control law identical to the one shown 
in the previous section. 

For this control law, the pseudocontrol, V 5 , is defined as 

^5 — hcom ( B-74) 

and the real nonlinear control is defined as in equation (B- 66 ) 

« = Pcom 

The intermediate real nonlinear control is defined as 

145 , = <f>com (B-75) 


The linearizing transformation used here is obtained from the vertical acceleration equation as a function of 
accelerometer outputs which is derived in equation (B-25) in the Derivation of Pitch Control Inverse Transformation 
Equations section 

h = a x sin 0 - a y sin <p cos 0 + o„ cos <f> cos 0 - 1 

This measurement equation is converted to a command equation by substituting 0 = <f> com and h = h com into it 
to get 

hcom = a x sin 0 - a y sin <j> com cos 0 + a n cos <f> com cos 0 - 1 ( B-76) 

From this and equations (B-75) and (B-74) the linearizing transformation is obtained as 


V 5 = a x sin 0 - a y sin us, cos 0 + a n cos us, cos 0 - 1 


(B-77) 
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The inverse linearizing transformation is obtained by solving equation (B-76) for 4> com using the trigonometric 
identity that given 

o cos 4> + b sin 4> = c 


with 

then 

which gives 


c 2 < o 2 + b 2 


, b _i c 

<f> = tan -1 - + cos 7 = ■ = 
a Vo 2 + b 2 


ficom — "*■ hin 


-1 


—a. 


0>t/ \ — 1 

y ’ + COS 


, + Ox sin 6 - 1 


(B-78) 

(B-79) 

(B-80) 

(B-81) 


cos 0 yoy 2 + a n 2 

Care mus. be taken with this equation to ensure that the angles obtained from eaeh of the terms fail in the correct 
quadrants. Using this equation and substituting v 5 = h com and u 5 , = An, '"1° 11 S'ves the inverse hncanzi g 
transformation 


us — — tan 


-l 


—a. 


J,y ' - cos 1 


-vs + a x sin 9 — 1 


cos 0 


•^Oy 2 + °" 2 


( B-82) 


The form of the pseudocontrol is identical to the form used in the pitch control of altitude shown in the Derivation 
of Pitch Control Inverse Transformation Equations section in equation (B-34) 


v = Kh = f kh, kh P kh D 


f eh 
eh 

£h 


(B-83) 


At this point, the commanded roll attitude, 4> CO m - is selected as 4> Te f ar| d 

= 4*00111 


( B-84) 


The inner layer control can then be determined using the roll attitude controller described in the Roll Attitude Com 
mand section. 

During gain selection, the inner layer of the outer-loop control can be similar, if not the same as those in the roll 
attitude controller. Again, the gains for the outer layer must be selected to produce slower dynamics than the inner 
layer to avoid invalidating the assumption of time-scale separation. 


Mach Number Command 

The Much number commanded roll control was not flight tested, however it was implemented and Med I in the 
simulation. It is used to track Mach number with roll attitude during maneuvers with fixed throttle at load factors 
greater than one, such as a thrust-limited windup turn. 

This control law is also a multilayer controller. In this case, however, the control law has lhrce la ^2' n .^ 
outermost layer generates an altitude rate command, then feeds it into an altitude commanded roll controller identical 
to the one discussed in the previous section which accounts for the other two layers in this controller. 

For this control law, the pseudocontrol is defined to be 


V6 = Me, 


(B-85) 
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while the real nonlinear control is again shown in equation (B-66) 


U ~ Pcom 

The first intermediate real nonlinear control is defined to be 

U6, = h com (B-86) 

The linearizing transformation is derived from equations (B-4) through (B-16). First, the velocity rate terms, 
ti, v, and w, are converted from functions of forces to functions of kinematic accelerations by substituting in equa- 
tions (B-6) through (B-8) to get 


u = 9a x ,k + rV sin /? - qV sin a cos (3 
v = ga y ,k + pV sin a cos j3 - rV cos a cos /? 
w = ga lik + qV cos a cos /? - pV sin /? 


(B-88) 

(B-89) 


Then, substitutions of accelerometer outputs for kinematic accelerations are made using equations (B-9) (B-10) 
(B-ll), and (B-12) to get 


** = 9( Q i — sin 9) + rV sin (3 — qV sin a cos (3 
v - 9( a y + sin <j) cos 9) + pV sin a cos f3 — rV cos a cos (3 
w = —g(a n - cos 0 cos 9) + qV cos acos^ - pVsin p 

These equations are substituted into equation (B-4) and terms are cancelled to get 

V = g(a x cos a cos [3 + a y sin/? - a n sin acos/? + 

sin </> cos 9 sin /? — sin 9 cos a cos /? + cos <f> cos 9 sin a cos /?) 

Next, equation (B-5) is solved for £ and is substituted into the gravity term of equation (B-93) giving 

V = 9 cos a cos /? + Oy sin (3 — o n sin a cos /3 — ^ j 
The derivative of equation (B-13) assuming constant acoustic velocity, is taken so that 

a 

and incorporated into equation (B-94) to give the measurement equation for M 

M = ~ fa* cos a cos /? + a y sin (3 - a„sin a cos & - 


(B-90) 

(B-91) 

(B-92) 


(B-93) 


(B-94) 


(B-95) 


(B-96) 


This measurement equation is converted to a command equation by setting M = M com and h = h com so that 

M com = £ cos a cos /? + a y sin {3 - a„ sin o-cos (3 - 

This equation is converted to the linearizing transformation by using equations (B-85) and (B-86) to get 


(B-97) 


V6 = a ( a * COSaCOS / 3 + a v sin (3 — a n sin a cos (3 — 


(B-98) 
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To get the inverse linearizing transformation, this equation is solved for U6, 

-v eaV 


U6, = 


• + V (a x cos a cos (3+ o y sin /3- o n sin a cos (3) 


(B-99) 


The pseudocontrol for this layer of the controller is simply a proportional-integral controller giving the follow- 
ingform r 

Jew | (B-100) 


V4 = [ Km,Km p ] 




The gains selected for this linear control law must not violate the separation of time scales for the inner layers. This 
control is. in general, a little bit slow and has difficulty maintaining precise control of Mach number. 

Derivation of Thrust Control Inverse Transformation Equations 

There were three controllers developed for thrust control: a Mach number controller, a velocity controller, and 
a dinmuhrodle »mmand. The direct tittle command is trivia, and is no, addressed here. The veloctty command 
is derived here, but was not tested in the simulation or in flight. 

The command required onboard the aircraft is PLA 

p r A ( B-101) 

U = PLAcom 

The inverse linearizing transforms used for this controller give an intermediate command value of thrust 

v ( B-102) 

tii = A Tcom 

which is converted to PLA C om using an inverse thrust model 

PI ,A cm = g(.X^,h,M) (B ' 103> 

The thrust model itself is used to estimate the thrust generated by the engines at the current PLAs 
is a simple version of the one used in the F-15 simulation. In this model. = f(PLA, , h 

duee PLA break points, three Mach number break points, and two alutude break points. The model uses linear 
interpolation between break points, and incorporates a simple lag, jfj, to model engine dynamics. 


Mach Number Command 

For this controller, the real nonlinear control is shown in equation (B-101) and the linear control is defined as 

VI - Mcom (B - 104) 

linearizing transform is derived from equation (B-4). To accomplish this the equations 1 for r ^elerometer 
of forces are derived by combining equations (B-6) through (B-8) with equations ( ) g 


The 

outputs as functions 
(B-12) which give 


a x = — ( Xt - D cos a + L sin a) 

gm 

a y =—Y 

gm 

a n = — ( -Zj + D sin a + L cos a) 

gm 


(B-105) 

(B-106) 

(B-107) 
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Next, equation (B-105) is multiplied by cos a, and equation (B-107) is multiplied by - sin a , and the two resulting 
equations are added together and reordered to get 

Zj sin a — D = gm(a x cos a — a„ sin a) — Xt cos a (B-108) 

For this equation, the accelerometer variables are measured, and the X -thrust variable is estimated with the thrust 
model so that Xy = Xr eat which gives 

Z T sma~D = gm( a x cos a — a n sin ck) — Xt c „ cos a ( B-109) 

This equation is substituted into equation (B-4) to eliminate the Zi sin a - D expression and equation (B-106) is 
substituted in to eliminate the Y -axis force terms giving 

• 1 

y ~ m 1 cos ™ 3 m ( °x cos a - 0. n sin a) - X Tej( cos a] + gma y sin 0+ X T cos a cos (3 + 


mg( cos ex cos (3 sin 9 — sin /?sin 0cos 9 — sin acos 0 cos 0 cos 9)} (B-110) 

The next step is to change the equation from V to M using equation (B-95) to get 
• _ 1 

~ cos 0(9 m ( a x cos a - a n sin a) - X Tej( cos a) + gma v sin 0 + X T cos acos (3 + 

— mg( cos a cos /? sin 9 — sin /3 sin cos 9 — sin a cos /? cos 0 cos 0) ] (B- 1 1 1 ) 

TOs equation is transformed from a measurement equation to a command equation by setting M = M com and 
= x r com to get 

■ _ 1 

com ~ am* C0S ^ gm ^ Qx C0S a ~ a " sin a - Xr '« cos + 9 ma y sin 0 + Xj com cos a cos (3 + 

— mg ( cos acos 0 sin 9 — sin /3 sin 0 cos 9 — sin acos /?cos 0cos 9) ] (B-l 12) 

By using equations (B-102) and (B- 104) this equation becomes the linearizing transform equation 
1 

Vl ~ am* C ° S cos a ~ a n sin a) — Xt c „ cos a) + gma y sin (3 + u, cos acos 0 + 

— mg ( cos acos sin 9 — sin /3sin 0cos 9 — sin acos /3cos 0 cos 0)] (B-l 13) 

which is subsequently solved for u, to get the inverse linearizing transform 

U ‘ = cos acos -p{ Viam - cos 01 9m(a x cos a - a n sin a) - X T „ t cos a] - gm 

( a y s > n 0 — cos acos ^sin 9 + sin /3sin 0cos 9 + sin acos /?cos 0cos 9)} (B-114) 

nonlinear contro1 1S obtained using the intermediate command value of thrust, u„ and the inverse engine 
model. This can be represented mathematically by combining equations (B-101), (B-102), and (B-l 03) to get 

u = g(ui,h,M) (B-l 15) 

The dynamics of the engine were not modeled in the inverse thrust model. 

Since this command is expected to be either a constant or a ramp, the linear control required is a proportional- 
integral control law similar to the one shown in equation (B-100) 


v 4 ~ [ KmjKmp J 


f cm 
eu 


The gains selected for this controller need to be slow enough to account for time-scale separation between the linear 
controller and the engine dynamics or the thrust model lag (they should be approximately the same). 
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Velocity Command 

For velocity control of PLA, equation (B-110) is used and the linear control is defined as 

Vg — Vcom 

The real nonlinear control remains the same as in equation (B-101) 

U = P LAcom 

and the intermediate control also remains the same 

U i = -^Tcon, 


(B-116) 


The linearizing transform, which is obtained from equation (B-l 10), would then be 


„ g = 1[ cos 0 (gm( ox cos a - a n sin a) - X TeJt cos a) + gma y sin 0 + u, cos a cos 0 + 
m 

-mg( cos a cos/? sin 0- sin 0 sin <j> cos 9 - sin a cos 0 cos <p cos 6) } 
and the inverse linearizing transform would be 


(B-l 17) 


(ug m - cos /?[ gm( a x cos a - o n sin a) - X Tcal cos a] - gm 

cos a cos (3 

( a y sin p - cos a cos 0 sin 9 + sin 0 sin <fr cos 6 + sin a cos 0 cos cf> cos 6) } 


(B-l 18) 


Since the velocity command is similar to the Mach number command, the control required is a proportional 
integral control law similar to the one shown in equation (B-100), except in terms of velocity so that 


The same limitations on gain 
this controller. 


V 4 


[ K v ,K Vp } 


f ey 
e v 


(B-119) 


selection as in the Mach number command in the previous section apply to 


49 


REFERENCES 


1. Burcham, Frank W., Jr., Gary A. Trippensee, David F. Fisher, and Terrill W. Putnam, Summary of Results of 
NASA F-15 Flight Research Program, NASA TM-8681 1, 1986. 

2. Menon, P.K.A., M.E. Badgett, R.A. Walker, and E.L. Duke, “Nonlinear Flight Test Trajectory Controllers for 
Aircraft,” J. Guidance, Control, and Dynamics, vol. 10, no. 1, Jan.-Feb. 1987, pp. 67-72. 

3. Menon, P.K.A., M.E. Badgett, and R.A. Walker, Nonlinear Maneuver Autopilot For The F-15 Aircraft NASA 
CR- 179442, 1985. 

4. Meyer, George and Luigi Cicolani, A Formal Structure For Advanced Automatic Flight-Control Systems 
NASA TN D-7940, 1975. 

5. Meyer, George and Luigi Cicolani, “Application of Nonlinear Systems Inverses to Automatic Flight Control 
Design System Concepts and Flight Evaluations,” AGARDograph AG-251 on Theory and Applications of 
Optimal Control in Aerospace Systems, edited by Pieter Kant, NATO, 1980, pp. 10-1 to 10-29. 

6. Meyer, George, “The Design of Exact Nonlinear Model Followers,” 1981 Joint Automatic Control Conference, 
vol. 2, Charlottesville, VA, June 17-19, 1981 , article FA-3A. 

7. Meyer, G„ R. Su, and L.R. Hunt, “Applications of Nonlinear Transformations to Automatic Flight Control,” 
International Federation of Automatic Control, Automatica, vol. 20, no. 1, 1984, pp. 103-107. 

8. Duke, Eugene L., Automated Flight Test Maneuvers: The Development of a New Technique,” 1982 Society 
of Flight Test Engineers (SFTE) Symposium, New York, NY, Sept. 20-22, 1982, pp. 101-119. 

9. Duke, Eugene L„ Michael R. Swann, Einar K. Enevoldson, and Thomas D. Wolf, “Experience with Flight 
Test Trajectory Guidance,” J. Guidance, Control, and Dynamics, vol. 6, no. 5, Sept.-Oct. 1983, pp. 393-398. 

10. Duke, E.L. and F.P. Jones, “Computer Control for Automated Flight Test Maneuvering,” J. Aircraft vol 21 
no. 10, Oct. 1984, pp. 776-782. 

11. Duke, Eugene L., Frank P. Jones, and Ralph B. Roncoli, Development and Flight Test of an Experimental 
Maneuver Autopilot for a Highly Maneuverable Aircraft, NASA TP-2618, 1986. 

12. Menon, P.K.A., R.A. Walker, and E.L. Duke, “Flight Test Maneuver Modeling and Control,” J. Guidance, 
Control, and Dynamics, A1AA-86-0426, Jan. 1986. 

13. Menon, P.K.A. and R.A. Walker, Aircraft Flight Test Trajectory Control, NASA CR-4161, 1988. 

14. Walker, Robert and Naren Gupta, Flight Test Trajectory Control Analysis, NASA CR-170395, 1983. 

15. Garg, Sanjay and David K. Schmidt, Optimal Cooperative Control Synthesis of Active Displays NASA 
CR-4058, 1987. 

16. Alag, Gurbux S. and Eugene L. Duke, “Development of Control Laws for a Flight Test Maneuver Autopilot,” 

J. Guidance, Control, and Dynamics, vol. 9, no. 4, July-Aug. 1986, pp. 441-^445. 

17. Menon, P.K.A., H.A. Saberi, R.A. Walker, and E.L. Duke, “Flight Test Trajectory Controller Synthesis with 
Constrained Eigenstructure Assignment,” 1985 4th American Control Conference, vol 3 Boston MA 
June 19-21, 1985. 

18. Antoniewicz, Robert F„ Eugene L. Duke, and Brian P. Patterson, User’s Manual for Interactive LINEAR, a 
FORTRAN Program to Derive Linear Aircraft Models, NASA TP-2835, 1988. 


50 



19. Duke, Eugene L„ Brian P. Patterson, and Robert F. Antoniewicz, User’s Manual for LINEAR, a FORTRAN 
Program to Derive Linear Aircraft Models, NASA TP-2768, 1987. 

20. Smith, G.A. and George Meyer, “Total Aircraft Flight-Control System — Balanced Open- and Closed-Loop 
Control with Dynamic Trim Maps,” IEEE Third Digital Avionics Systems Conference, Fort Worth, TX, Nov. 

6-8, 1979, pp. 215-223. 

2 1 . Smith, G. Allan and George Meyer, Application of the Concept of Dynamic Trim Control to Automatic Landing 
of Carrier Aircraft, NASA TP-1512, 1980. 

22. Smith, G. Allan and George Meyer, “Application of the Concept of Dynamic Trim Control and Nonlinear Sys- 
tem Inverses to Automatic Control of a Vertical Attitude Takeoff and Landing Aircraft,” 
AIAA-8 1-2238, 1981. 

23. Smith, G. Allan and George Meyer, “Aircraft Automatic Digital Flight Control System With Inversion of the 
Model in the Feed-Forward Path,” AIAA-84-2627, 1984. 

24. Smith, G. Allan and George Meyer, “Aircraft Automatic Flight Control with Model Inversion, J. Guidance, 
Control, and Dynamics, vol. 10, no. 3, May-June 1987, pp. 269-275. 

25. Wehrend, William R., Jr., Pilot Control through the TAFCOS Automatic Flight Control System, NASA 
TM-81152, 1979. 

26. Wehrend, William R., Jr. and George Meyer, Flight Tests of the Total Automatic Flight Control System 
(TAFCOS) Concept on a DHC-6 Twin Otter Aircraft, NASA TP-1513, 1980. 

27. Meyer, G., R.L. Hunt, and R. Su, Design of a Helicopter Autopilot by Means of Linearizing Transformations, 
NASA TM-84295, 1982. 

28. Hunt, L.R., R. Su, and G. Meyer, “Approximating Linearizations for Nonlinear Systems, NASA TM-88772, 
1986. 

29. Hunt, L.R. and Renjeng Su, “Approximations of Nonlinear Systems Having Outputs,” NASA CR- 176946, 
1984. 

30. Hunt, L.R., Renjeng Su, and George Meyer, “Global Transformations of Nonlinear Systems,’ IEEE Transac- 
tions on Automatic Control, vol. AC-28, no. 1, Jan. 1983, pp. 24-31. 

31. Meyer, George, Renjeng Su, and L.R. Hunt, Applications to Aeronautics of the Theory of Transformations of 
Nonlinear Systems, NASA TM-84249, 1982. 

32. Su, Renjeng, L.R. Hunt, and George Meyer, Canonical Forms for Nonlinear Systems, NASA 
TM-81323, 1983. 

33. Su, Renjeng and L.R. Hunt, A Canonical Form for Nonlinear Systems, NASA CR-176974, 1986. 

34. Duke, Eugene L„ Robert F. Antoniewicz, and Keith D. Krambeer, Derivation and Definition of a Linear 
Aircraft Model, NASA RP-1207, 1988. 

35. Kato, Osamu and Ichiro Sugiura, “An Interpretation of Airplane General Motion and Control as Inverse Prob- 
lem,” J. of Guidance, Dynamics, and Control, vol. 9, no. 2, Mar. -Apr. 1986, pp. 198-204. 

36. Cicolani, Luigi S. and Stein Weissenberger, A Nonlinear Trajectory Command Generator for a Digital Flight- 
Control System, NASA TP- 1221, 1978. 


51 


REPORT DOCUMENTATION PAGE 


1. AGENCY USE ONLY (Lnv blink) 


TITLE AND SUBTITLE 


REPORT DATE 

September 1991 


Form Approved 
OMB No. 0704-0188 


ravwwing Instruction*. Marching ailstlng data lourcn, 
regarding this burdan animat# or any othar aapact of this 
la for inform* ion Operation* and Raporta. 1 215 Jaftaraon 
>n Profact (0704-01 SI). Waahlngton. OC 20503. 


REPORT TYPE ANO DATES COVERED 
NASA Technical Paper 

js. FUNDING NUMBERS 


Application and Flight Test of Linearizing Transformations Using 
Measurement Feedback to the Nonlinear Control Problem 


6. AUTHOR(S) 

Robert F. Antoniewicz, Eugene L. Duke, and P.K.A. Menon 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Dryden Flight Research Facility 
RO. Box 273 

Edwards, California 93523-0273 


fl. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


11. SUPPLEMENTARY NOTES 1 

Robert F. Antoniewicz and Eugene L. Duke: Dryden Flight Research Facility, 
Edwards, California. 

P- K- A. Menon: Georgia Institute of Technology, Atlanta, Georgia, 


RTOP 505-60-21 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


H-1629 


1 10. SPONSORINQ/MONITORING 
AGENCY REPORT NUMBER 


NASA TP-3154 


■ 1*«. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified — Unlimited 
Subject Category 08 

13 - ABSlRAC^fUMxhmjrrT^OO 1 , m- 

The design of nonlinear controllers has relied on the use of detailed aerodynamic and engine models that must 
be associated with the control law in the flight system implementation. Many of these controllers have been 
applied to vehicle flightpath control problems and have attempted to combine both inner- and outer-loop control 
functions in a single controller. In this paper, a new approach to the nonlinear trajectory control problem is 
presented. This approach uses linearizing transformations with measurement feedback to eliminate the need for 
detailed aircraft models in outer-loop control applications. By applying this approach and separating the inner- 
loop and outer-loop functions two things are achieved: (1) the need for incorporating detailed aerodynamic 
models in the controller is obviated and (2) the controller is more easily incorporated into existing aircraft flight- 
control systems. This paper also discusses an implementation of the controller described here. This controller 
is tested on a six degree -of-fteedom F-15 simulation and in flight on an F-15 aircraft Simulation data are 
presented which validates this approach over a large portion of the F-15 flight envelope. Proof of this concept 
is provided by flight-test data that closely matches simulation results. Right-test data are also presented. 


14. SUBJECT TERMS "" 

15. NUMBER OF PAGES 

Aircraft control; Dynamic conversions; Linearizing transformations; 

56 

Nonlinear control 

16. PRICE CODE 

A04 


OF THIS PAGE 

Unclassified 


OF ABSTRACT 

Unclassified 


Standard Form 298 (Rev. 2-09) 

PrM«rK>*d by ANSI Std. Z34-1 8 
2*4-1 02 


NAS A -Langley. 1991 







National Aeronautics and 
Space Administration 
Code NTT-4 


Washington, O.C. 




01 he rat Business . 

; Penalty (Or Private Use, $300 



BUCK RATE 

POSTAGE & FEES PAID 

NASA 

Permit No. G-27 , - 


POSTMASTER: 


If UixS«Hv«nifel+ {Suction 158 
Postil Manual) Do Not Return 



